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Abstract 

Trypanosoma cruzi, the causative agent of Chagas disease, is extremely resistant to ionizing radiation, enduring up to 
1.5 l<Gy of gamma rays. Ionizing radiation can damage the DNA molecule both directly, resulting in double-strand breaks, 
and indirectly, as a consequence of reactive oxygen species production. After a dose of 500 Gy of gamma rays, the parasite 
genome is fragmented, but the chromosomal bands are restored within 48 hours. Under such conditions, cell growth 
arrests for up to 120 hours and the parasites resume normal growth after this period. To better understand the parasite 
response to ionizing radiation, we analyzed the proteome of irradiated (4, 24, and 96 hours after irradiation) and non- 
irradiated r. cruzi using two-dimensional differential gel electrophoresis followed by mass spectrometry for protein 
identification. A total of 543 spots were found to be differentially expressed, from which 215 were identified. These 
identified protein spots represent different isoforms of only 53 proteins. We observed a tendency for overexpression of 
proteins with molecular weights below predicted, indicating that these may be processed, yielding shorter polypeptides. 
The presence of shorter protein isoforms after irradiation suggests the occurrence of post-translational modifications and/or 
processing in response to gamma radiation stress. Our results also indicate that active translation is essential for the 
recovery of parasites from ionizing radiation damage. This study therefore reveals the peculiar response of 7. cruzi to 
ionizing radiation, raising questions about how this organism can change its protein expression to survive such a harmful 
stress. 
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Introduction 

Chagas disease, a neglected tropical disease caused by the 
protozoan parasite Trypanosoma cruzi, is considered to be a public 
health problem [1,2]. Over 10 million people are infected in Latin 
America and more than 100 million individuals live at risk of 
infection by blood transfusion, congenital, or oral transmission [3] . 
Forty years after its introduction, benznidazole and nifurtimox 
continue to be the first choice of treatment for Chagas disease. 
However, chemotherapy based on nitroheterocyclic compounds 
has a limited efficacy for patients in the chronic phase of infection 
and these drugs are highly toxic [4,5]. Little progress has been 
made toward the treatment of infected individuals and the 
development of more efficient drugs to treat Chagas disease 
patients remains urgent. Considering the resistance of some 
parasites to chemotherapy, the introduction of vaccines against T. 
cruzi could be another option [3,6] . 

T. cruzi is capable of resisting high doses of gamma radiation, 
enduring up to 1.5 kGy. As a direct biological effect, gamma 
radiation causes double-strand breaks (DSB) in the parasite DNA. 
However, 48 hours after irradiation, it is possible to see the 
chromosomal bands already restored. The parasite growth arrests 



for up to 120 hours, returning to the normal rate after this period 
[7,8] . This extraordinary recovery might be due to a very efficient 
DNA repair system. Homologous recombination is required to 
repair DNA DSBs and the involvement of the TcRAD51 protein 
in this process was evaluated by our group elsewhere. The 
overexpression of TcRAD51 ensures a more effective DSB DNA 
repair and a greater resistance to DNA damage in T. cruzi [9] ■ 

Oxidative stress is another effect of ionizing radiation due to the 
production of hydroxyl radicals (OH'), superoxide (O2"), and 
hydrogen peroxide (H2O2), directly from radiolysis of water. These 
products are commonly called reactive oxygen species (ROS) [10]. 
Once the DNA molecule is intimately associated with water, the 
production of OH' results in damages that include, apart from 
DSBs, oxidation of nitrogenous bases and sugar [11,12]. 
Approximately 75-80% of the biological damage caused by this 
type of radiation is mediated by OH" formation. Such radicals are 
capable of reacting with most biologically relevant molecules. Each 
amino acid reacts differendy with OH' and the precise mecha- 
nisms of reaction are poorly understood [13]. 

Another organism that is extremely resistant to ionizing 
radiation is the bacterium Deinococcus radiodurans, which can 
withstand radiation doses of up to 15 kOy [14]. D. radiodurans 
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presents a very robust DNA repair apparatus; nevertheless, the 
biological responses to genomic lesions depend on its proteome 
integrity. Considering that ionizing radiation also induces protein 
damage through oxidative stress, a protected functional proteome 
ensures an efficient cell recovery from this type of stress [15]. 
Using the classical proteomic approach of two-dimensional 
differential gel electrophoresis (2D-DIGE) coupled with mass 
spectrometry (MS), Basu & Apte observed in a time-course 
analysis that some classes of proteins have a strong influence on 
stress responses. These proteins are mainly involved in processes 
such as DNA damage repair, protein synthesis and folding, and 
responses to oxidative stress [16]. 

Proteome versus transcriptome analyses have been highly 
recommended for studies with tripanosomatids, as they have very 
peculiar molecular features concerning their gene expression 
control. As a kinctoplastid, T. cruzi transcription is polycistronic 
and gene regulation occurs mainly post-transcriptionally, with 
mature mRNAs being generated by trans-splicing and polyade- 
nylation [17,18]. The processing and stabilization of mRNAs are 
extremely important in trypanosomatid gene regulation [19,20]. 
Furthermore, other dynamic control mechanisms, such as post- 
translational modifications, are fundamental in the regulation of 
gene expression and need to be better characterized in these 
organisms [21-23]. 

A time-course microarray study previously carried out by our 
group analyzed the T. cruzi gene expression in response to gamma 
radiation [7]. Among the 273 dilferentially expressed genes, 160 
were upregulated and 1 1 3 were downregulated. The majority of 
the genes with assigned functions was downregulated. Translation, 
protein metabolic processes, and the generation of precursor 
metabolites and energy pathways were affected. Four mitochon- 
drial genes and Retrotransposon Hot Spot genes were upregu- 
lated; likewise, the tyrosyl-DNA phosphodiesterase 1, a gene 
involved in DNA DSB repair, was also induced [7]. Taking into 
account the T. cruzi gene expression peculiarities, analyses of 
proteome changes after irradiation in different time points may 
contribute to the understanding of the parasite response to such 
stress. 

In this work, we performed quantitative proteomic analyses 
using 2D-DIGE to ascertain the parasite response to ionizing 
irradiation. A total of 543 protein spots were found to be 
differentially expressed considering all analyzed time points and 53 
difiFerent proteins were identified by tandem mass spectrometry 
(MS/MS). The great majority of the identified proteins was 
represented by several isoforms, suggesting that post-transcrip- 
tional and/ or post-translational modifications are occurring as a 
consequence of gamma radiation exposure. Overexpression of 
tryparedoxin after irradiation was also observed, indicating that 
the parasite may be responding to the oxidative stress caused by 
irradiation. We also compared the time-course microarray and 
proteomic analyses. Although some of the protein expression 
patterns confirmed the microarray results, the correlation between 
mRNA and protein levels of the genes identified in both studies 
was extremely poor. In addition, treatment of the parasites with 
translation inhibitors showed that the synthesis of proteins 
putatively involved in the parasite response to stress is essential 
for its recovery from such a harmful stress. 

Materials and Methods 

Cell Culture and Gamma Irradiation 

In this work, we used T. cruzi epimastigote forms of the CL 
Brener strain, which were isolated and characterized by Brener & 
Chiari [24]. Clones have been maintained as frozen stocks at 



Universidade Federal de Minas Gerais. Parasites were grown at 
28°C in liver infusion tryptose (LIT) medium pH 7.3, supple- 
mented with 10% fetal bovine serum, streptomycin sulfate (0.2 g/ 
L), and penicillin (200,000 units/L). Cultures in the exponential 
growth phase (2x10 cells/mL) were exposed for 20 minutes to 
500 Gy of gamma radiation (1,578 Gy/h) in a cobalt (60 Co) 
irradiator (Centre de Desenvolvimento da Tecnologia Nuclear - 
CDTN, Belo Horizonte, Brazil). Cells were counted daily after 
irradiation to generate the growth curve. 

Cycloheximide and Puromycin Treatments 

Parasites exposed or not exposed to 500 Gy of gamma radiation 
were treated with cycloheximide (Calbiochem) 50 |ig/mL for 15 
minutes or with puromycin (Sigma) 25 |ig/mL for 1 hour. Both 
drugs were added to the parasite cultures 4 hours after irradiation. 
Parasites were washed twice in phosphate buffered saline (137 mM 
NaCl, 4 mM NagHPO*, 1.7 mM KH2PO4, and 2.7 mM KCl), 
the LIT medium was replaced, and the cells were counted. 

Protein Extract Preparation and DIGE Labeling 

Protein extracts were obtained, simultaneously, in triplicate for 
each condition: non-irradiated control (NI), 4, 24, and 96 hours 
after irradiation. Parasites (2x10^ cells) were washed twice with 

LIT medium followed by centrifugation at 1,500 g for 5 minutes at 
4°C. Each pellet was resuspended in 200 |J.L of lysis buffer (8 M 
urea, 2 M thiourea, 4% CHAPS, 10 mM Tris base) and a 
protease inhibitor mix (GE Healthcare, USA). Samples were 
mixed on vortex every 30 minutes during 2 hours of incubation at 
room temperature and subsequentiy centrifuged at 14,000 g for 30 
minutes. The supernatants were aliquoted and stored at — 70°C 
for further use. For all samples, protein concentration was 
determined using the 2D Quant kit (GE Healthcare, USA), 
according to manufacturer's instructions. 

Before labeling, samples had their pH adjusted to 8.5 with 
NaOH 0.05 M (as recommended by the manufacturer's protocol). 
To reduce biological variation, a pool of protein (extracts obtained 
from triplicates was used. A total of 50 ng of protein from each 
pool (NI, 4, 24, and 96 hours after irradiation) was labeled with 
CyDye DIGE Fluor Minimal Labeling Kit (GE Healtiicare, USA). 
The dye swap strategy was used to avoid label bias, where each 
sample was labeled with 400 pmol of either Cy3 or Cy5. A 
mixture of all protein extracts (12.5 |J.g of each pool sample) was 
labeled with Cy2 as the internal control. Reactions were carried 
out on ice for 30 minutes in the dark and then stopped by the 
addition of 10 mM lysine. 

Two-Dimensional Gel Electrophoresis 

First dimension. The isoelectric focusing (lEF) was per- 
formed using Immobiline Dry Strips (GE Healthcare, USA) 1 8 cm 
in size, with a pH ranging from 4-7. Strips were loaded with 50 (ig 
of protein per CyDye (total of 150 |j.g) and sample buffer 
containing 8 M urea, 2 M thiourea, V/o CHAPS, 1% dithiothre- 
itol (DTT), 0.002% bromophenol blue, and \% IPG buffer 
(pH 4—7; GE Healthcare, USA). Passive rehydration followed 
overnight, at room temperature, in a strip holder (GE Healthcare, 
USA). The IFF protocol used in the Ettan IPGphor3 (GE 
Healthcare, USA) instrument was as follows: 50 |i,A per strip, 
20°C, steps 1 to 5: 0.2 kV for 12 hours, 0.5 kV for 2 hours; 1 kV 
for 1.5 hour, 8 kV for 2 hours, 8 kV gradually raising to 40 kV, 
accumulating approximately 60 kV in total. Focused IPG strips 
were equilibrated for 15 minutes in an equilibration solution 
(50 mM Tris-HCl pH 8.8, 6 M urea, 30% glycerol, 2% SDS, 
0.002% bromophenol blue and 125 mM DTT) and tiien alkylated 
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for an additional 15 minutes in an equilibration solution 
containing 13.5 mM iodoacetamide instead of DTT. 

Second dimension. Equilibrated strips were briefly washed 

in Ix running buffer (25 mM Tris, 192 mM glycine, and 0.2% 
SDS) and placed on top of 12% acrylamide/bis-acr^-lamide gels, 
overlaid with a 0.5% agarose solution. Protein separation was 
carried out at 10°C, in an Ettan Dalt Six Electrophoresis System 
(GE Healthcare, USA), 45 mA per gel, until the dye front reached 
the bottom of the gel. Labeled proteins in each gel were visualized 
using the Typhoon FLA 9000 scanner (GE Healthcare, USA) at 
100 |a,M image resolution with excitation/emission wavelengths 
for Cy3 (532/580 nm), Cy5 (633/670 nm), and Cy2 (488/ 
520 nm). Gel images were uploaded and cropped using Image 
Loader Software (GE Healthcare, USA), then imported to 
DeCyder 2D software, version 7.0 (GE Healthcare, USA). 

DIGE Data Analysis 

For spot detection, the Differential In-gel Analysis (DIA) module 
of DeCyder 2D software, version 7.0 (GE Healthcare, USA), was 
used. The DIA co-detection algorithm exploits the identical spot 
patterns from multiple samples in the same gel. After the removal 
of some artifacts from the gels, spot quantification was performed 
automatically by normalizing the spot volumes against the internal 
control. The following steps were performed in the Biological 
Variation Analysis module, which uses images processed in DIA 
and matches spots across gels. One-way ANOVA and Student's t- 
test were applied to evaluate differential protein expression levels 
between the groups of study. Spots classified as significandy 
differentially expressed were manually inspected. Abnormal spots 
were excluded from the analysis when necessary and gels were re- 
matched. 

Trypsin in-Gel Digestion, Mass Spectrometry, and Protein 
Identification 

DifiFerentially expressed protein spots were excised and trypsin 
in-gel digestion was carried out overnight at 37°C with 20 ng/|j,L 
of trypsin (Promega, Sequencing Grade Modified Trypsin, USA), 
diluted in 25 mM ammonium bicarbonate. After trypsin digestion, 
peptides were extracted from the gel by washing twice with 30 |J,L 
of 50% acetonitrik- and 5% formic acid solution and shaking for 
15 minutes. Peptides were then concentrated (Eppendorf Con- 
centrator 5301) to 10 [lL and desalted using Zip-Tip (CIS resin, 
PIO, MiUipore Corporation, USA). Once the peptides were eluted 
(50% acetonitrile/0.1% trifluoroacetic acid) from columns, 0.5 |J,L 
of each sample was mixed with 0.25 |J,L of a saturated matrix 
solution [10 mg/mL a-cyano-4-hydroxycinnamic acid (Aldrich, 
USA) in 50% acetonitrUe/O.l'/'o trifluoroacetic acid]. Samples 
were spotted on the MTP AnchorChip 600/384 (Bruker 
Daltonics) and let to dry at room temperature. Raw data for the 
identification of proteins were obtained with the MALDI-TOF- 
TOF AutoFlex III (Bruker Daltonics, USA) instrument (Labor- 
atorio Multiusario de Biomoleculas, Departamento de Bioquimica 
e Imunologia, UFMG, Brazil) in the positive/reflector mode 
controlled by FlexControl software. Instrument calibration was 
achieved by using peptide calibration standard II (Bruker 
Daltonics) as a reference. Trypsin and keratin contamination 
peaks were excluded from the peak lists used for data base 
searching. Each spectrum was produced by accumulating data 
from 200 consecutive laser shots. 

MS/MS spectra were searched against the non-redundant 
protein sequence database from the National Center for Biotech- 
nology Information (http://www.ncbi.nlm.nih.gov) using the 
MASCOT software (version 2.1) MS/MS ion search tool 
(http://www.matrixscience.com). The search parameters were as 



follows: no restrictions on protein molecular weight, two tryptic 
miss-cleavages allowed, and variable modifK:ations of methionine 
(oxidation), cysteine (carbamidomethylation), and pyroglutamate 
formation at N-terminal glutamine of peptides. The mass 
tolerance for the peptides in the searches was 0.6 Da for MS 
spectra and 0.4 Da for MS/MS spectra. Peptides were considered 
to be identified when the scoring value exceeded the identity or 
extensive homology threshold value calculated by the MASCOT 
software (p<0.05). 

Manual Curation and Statistical Analysis 

Peptide sequences obtained from MASCOT were ahgned to the 
T. cruzi annotated genome using the BLAST tool from TriTrypDB 
(http://www.tritrypdb.org). Protein annotation was reassigned 
particularly when partial sequences were chosen by MASCOT 
and full-length sequences were available at the TriTr)'pDB. Once 
a final annotated and curated set of upregulated and downreg- 
ulated spots was available, it was possible to assess the protein 
species by their expected and observed weights (retrieved from the 
TriTrypDB and calculated from the position in the 2D-DIGE, 
respectively). 

Statistical analyses were performed using R in-house scripts with 
built-in statistical functions. A linear model was apphed to test the 
correlation between molecular weight and fold-change. The 
Wilcoxon test was used to evaluate the presence of significant 
differences between 1) the observed molecular weights of 
upregulated and downregulated protein spots and 2) the observed 
and (-xpcctc-d mok-cular weights from upregulated and downreg- 
ulated protein spots. 

The final set of proteins was further manually annotated 
according to biological function and grouped into difierent 
functional classes based on literature data describing each protein 
and its molecular role. 

Results and Discussion 

The Effects of Protein Synthesis Inhibition on the Growth 
of T. cruzi Epimastigote Cells Exposed to Gamma 
Radiation 

Normal growth of epimatigote cells was affected by protein 
synthesis inhibition (using 50 Hg/ mL cycloheximide or 25 (Xg/mL 
puromycin) and by ionizing radiation treatment (500 Gy), as 
shown in Figure 1 . However, irradiation promoted a more drastic 
growth arrest that persisted for approximately 96 hours; after this 
period, the parasites resumed normal growth, reaching the 
stationary phase 216 to 240 hours after irradiation (Figure 1). 
The treatment of NI cells with cycloheximide (Figure lA) or 
puromycin (Figure IB) retarded the cell growth by at least 
24 hours when compared with non-treated cells, but did not lead 
to parasite death. Conversely, the combination of cycloheximide 
treatment and gamma radiation was lethal to 40% of the parasites. 
The remaining parasites resumed growth only 270 hours after 
irradiation, reaching the stationary phase 408 hours after irradi- 
ation (Figure lA). For puromycin, a similar effect was observed, 
but treated cells resumed normal growth eariier when compared 
with cycloheximide-treated parasites (Figure IB) and, in this case, 
no parasite death was detected. 

These observations indicate that an active translation is 
important for the recovery of parasites from damage caused by 
ionizing radiation. Protein synthesis blockage is potentially 
impairing the translation of newly synthesized or pre-existing 
mRNAs that code for proteins involved in triggering cell 
proliferation. These proteins may accumulate within 24 hours 
after irradiation and act later after irradiation, when parasites have 
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Figure 1. The effect of irradiation and translation inlilbltion on T. cri/z/ epimastigotes growth. Irradiated (500 Gy) or Nl parasites were 
treated with cycloheximide 50 |ig/mL (A) or puromycin 25 ng/mL (B), both added 4 hours after irradiation. Each point represents the mean ± 
standard deviation of three different experiments. 
doi:1 0.1 371/journal.pone.0097526.g001 



their genome recomposed and are ready to resume growth. 
Interestingly, our previous worij using microarray did not 
demonstrate upregulation of a variety of mRNAs coding for 
proteins involved in cell proliferation or cell cycle control, even 
96 hours after irradiation, and only a cyclin transcript was found 
to be affected (1.67-fold change increase) 24 hours after irradiation 
[7]. 

Analysis of the Proteome Profile of T. cruzi Epimastigote 
Cells Exposed to Gamma Radiation 

Since we have verified that newly synthesized proteins lia\ (; an 
impact on parasite recovery from irradiation stress, we decided to 
analyze time-course T. cruzi changes in the proteome induced by 
irradiation. Protein extracts were obtained from control NI cells 
and 4, 24, and 96 hours after irradiation. No significant losses in 
the total protein content and integrity were observed by ID-gel 
electrophoresis (Figure SI). Using the 2D-DIGE approach, six gels 
were produced following the experimental design specified in 
Table 1. This technique was chosen due to its greater sensitivity. 



reduced gel-to-gel variation, and its capacity for quantitative 
measurements of the relative abundance of each protein in a 
complex sample [25] . Figure 2 illustrates 2D-DIGE gels at all time 
points. An average of 2,1 86 ± 140 spots was found when compared 
with the master gel. From those, 543 presented altered expressions 
after irradiation, considering all time points (one-way ANOVA, 
p<0.01) and 215 were identified by peptide mass fingerprint, 
corresponding to 53 different proteins (Table 2). Almost half of 
these proteins (26) were represented by more than one spot in the 
2D gel (ranging from 2-12 spots per protein), indicating the 
presence of several isoforms for the same protein. These results 
suggest that post-translational modifications or protein processing 
are occurring during the response to gamma radiation stress. We 
have manually annotated the function of all 53 identified proteins 
via a literature search. Proteins were then manually assigned to 1 5 
different classes according to their biological function (Figure S3). 

Additionally, the Student's t-test (p<0.01) was applied to verify 
which proteins were differentially expressed in each time point 
when compared with the NI sample. The overall and time-specific 
number of downregulated protein spots was higher than the 
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Table 1. Experimental design. 
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Each two-dimensional gel was loaded with 50 ng of total protein extract per 
sample, labeled either with Cy3 or Cy5. The internal control (a pool containing 
50 Jig of all time point proteins: Nl, 4, 24, and 96 hours after irradiation) was 
labeled with Cy2. 

doi:l 0.1 371 /journal.pone.0097526.t001 



number of upregulated ones (Figure 3A). These findings are 
different from those described previously in our microarray study. 
Twenty-four hours post-irradiation, the number of downregxilated 
genes decreases drastically, reaching only 6 down-expressed 



transcripts 96 hours after irradiation, while the number of 
upregulated genes increases [7]. A linear regression analysis 
between mRNA and protein levels from genes concomitantiy 
identified in both studies was carried out for each time point. The 
correlation was extremely poor at all time points, starting with 
multiple R2 = 0.064 at 4 hours and reaching R2 = 0.27 and 0.24 
at 24 and 96 hours post-irradiation, respectively (data not shown). 
Although a very low correlation was obtained, the result is in 
agreement with other studies performed in both prokaryotes and 
eukaryotes using classical methodologies such as microarray. Serial 
Analysis of Gene Expression (SAGE) and RNA-Seq for transcrip- 
tomic expression data, and 2-DE, Multi-dimensional protein 
identification technology, and MS for proteomics data [26-31]. 
This analysis reinforces the idea that transcriptomic and proteomic 
approaches are complementary, not confirmatory [29] . 

Moreover, changes in the T. cruzi proteome are more evident 
24 hours after exposure to gamma radiation (Figure 3). This 
scenario suggests that epimastigote cells present an immediate but 
subde response to gamma radiation characterized by 12 induced 
and 21 repressed protein spots 4 hours after irradiation 
(Figure 3A). Between 4 and 24 hours after irradiation, a more 
intense response to stress was observed and most of the induced 




Figure 2. 2D-DIGE analysis of total protein extracts of irradiated and Nl epimastigote cells. Gel images 1-6 (see the experimental design 
in Table 1) showing - in triplicate - parasite proteins from each time point, labeled either with Cy3 (green) or Cy5 (red). Proteins were separated in the 
first dimension along a pH gradient (pH 4-7, 18 cm Immobiline DryStrip (GE Healthcare, USA), and in the second dimension in a 12% polyacrylamide 
gel. The molecular weight marker (MW) is indicated in kDa. 
doi:1 0.1 371/journal.pone.0097526.g002 
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Hours post-irradiation 



Figure 3. Protein spots differentially expressed at all time points. A) Number of downregulated and upregulated protein spots per time 
point. B) Venn diagram showing the overlaps of 32 protein spots differentially regulated among the three time points and of the 428 protein spots 
between 24 and 96 hours. 
dol:10.1371/journal.pone.0097526.g003 



and repressed protein spots were still significantly altered until 
96 hours (428 spots; Figure 3B). This finding indicates a sustained 
alteration in the abundance of specific T. cmzi proteins 24 hours 
after gamma radiation exposure. When analyzing the T. cmzi 
proteome 24 hours after irradiation, we found that, from the 59 
exclusive spots, approximately 66% were repressed and 34% were 
induced. However, the majority of the 23 exclusive spots found 
96 hours after irradiation were induced (approximately 61%). 

Exposure to Gamma Rays Increases the Levels of Shorter 
and/or Processed Proteins in Epimastigote Cells 

When analyzing the set of upregulated proteins (especially 24 
and 96 hours after irradiation), we observed a tendency for the 
overexpression of shorter molecules to the detriment of longer 
ones. The upregulated protein spots (red-colored dots) are mainly 
at the lower part of the gel (lower molecular weight), while the 
downregulated protein spots (green-colored dots) are more 
sparsely distributed across the gel (Figure 4A). In addition, low 
molecular weight protein spots tended to have larger fold changes 
when compared with those with molecular weights close to the 
expected value (Figure 4B). The WUcoxon test was applied and 
confirmed that the median values of the molecular weight of 
downregulated and upregulated protein spots were different for 
each time point (p<le-09; median values of 55.45/19.39, 45.64/ 
19.38, and 46.51/19.42 for 4, 24, and 96 hours, respectively; 
Figure 5A). 

When we considered the expected molecular weight of the fuU- 
length isoform (predicted size obtained from the TriTryDB 
website) of both upregulated and downregulated proteins, we 
noticed a decrease in protein size in the former case and an 
increase in the latter case, thus showing the emergence of lower 
molecular weight protein isoforms after irradiation (Figure 5B). 
We then decided to confirm if the observed molecular weight of 
these proteins in the 2D gels was in agreement with their expected 
molecular weight. In the case of upregulated proteins, the 
observed molecular weight was significandy lower than expected 
(Figure 5B), indicating that these proteins might be processed, 
yielding shorter polypeptides. It is important to note that this result 
does not seem to be a consequence of protein degradation, since 
clear spots can be observed in the 2D gel, indicating the presence 
of a large amount of identical polypeptides in this region of the gel. 



This would not be the case if proteins were degraded, considering 
that in this situation peptides of variable size would be generated 
and no clear spot would be observed in the gel. These results may 
indicate the emergence of new protein isoforms, as the result of 
protein processing, alternative splicing of mRNAs, and/or 
alternative translational start/stop sites after irradiation. Alterna- 
tive splicing of transcripts has the potential to expand the 
repertoire of proteins. Recent studies have estimated that all 
multi-exonic human genes are able to produce at least two 
alternatively spliced mRNA transcripts by alternative splicing, 
generating different proteins isoforms with altered structures and 
biological functions [32] . In trypanosomatids, mature mRNAs are 
generated after two processing events: trans-sphcing to add the 
spliced leader (SL) sequence to the 5' end of transcripts and 
subsequent polyadenylation [17]. A genome-wide analysis com- 
paring the SL addition site along the developmental cycle of the 
parasite suggests that alternative trans-sphcing plays an important 
role in differential gene expression [33]. The occurrence of 
alternative trans-splicing could be an explanation for the presence 
of so many different isoforms in T. cruzi after radiation response. 

A similar event has already been described in D. mdiodumns, 
since different isoforms of the single-strand binding (SSB) protein 
were produced after ionizing radiation stress induction. SSB 
proteins are vital for cell survival due to their involvement in 
processes such as DNA replication, recombination, and repair. 
The SSB protein spots in the gel followed a dynamic pattern of 
appearance, indicating a progressive processing of the C-terminal 
acidic tail, perhaps upon its interaction with ssDNA. The observed 
isoelectric point (pi) and molecular weight of deinococcal SSB 
isoforms were in agreement with the in jz&o-predicted pi and 
molecular weight of the SSB proteins shortened from the C- 
terminal end [16]. 

Intriguingly, in most of the observed processed proteins, the 
identified peptide sequences were the same or nearly the same in 
all sequenced protein spots and, therefore, it was impossible to 
define the actual outcome of the protein processing. As a 
particular case of study, the protein annotated as prostaglandin 
F2 alpha synthase, which is similar to NADH-flavinoxidoreduc- 
tase, is processed to a total of six different forms (Figure S2). While 
the expected molecular weight for the annotated sequence is of 
42 kDa, only two isoforms are nearly this size (~40 kDa) and 
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Figure 4. Distribution of upregulated and downregulated protein spots versus molecular weight, pi, and fold change. In the scatter 
plots, upregulated protein spots are shown in red and downregulated protein spots are shown in green. The correlation between molecular weight 
and pi or fold-change ratio is shown in (A) and (B), respectively. Spots with no significant difference in expression are colored gray. The blue line 
indicates the negative correlation between molecular weight and fold change. 
doi:1 0.1 371 /journal.pone.0097526.g004 



were, in fact, the most downregulated isoforms. A third isoform 
has a predicted molecular mass of 68 kDa, greatly exceeding the 
expected protein size. As all of the MS/MS-identified peptides 
were mapped to the C-terminal portion of the isoforms, there is no 
information to characterize the N-terminal of this enlarged protein 
naturally present in the NI parasites and downregulated after 
exposure to gamma radiation. A smaller (29 kDa) protein is 
expressed in approximately equal levels before and after radiation 
exposure, while an even smaller (22 kDa) protein species is 
exclusively present in irradiated cells. This is an interesting 
example, representative of multiple cases, in which we have 
observed the emergence of shorter isoforms of a same protein after 
epimastigote irradiation. The list of processed proteins expressing 
shorter isoforms after irradiation includes alpha and beta-tubulin, 
D-isomer specific 2-hydroxyacid dehydrogenase-protein, elonga- 
tion factor 2, glycerate kinase, pyruvate dehydrogenase El beta 
subunit, tyrosine aminotransferase, and several heat shock proteins 
(HSPs), such as HSP60, DnaK, HSP70s, and glucose-regulated 
protein 78. Apart from the previously discussed SSB proteins in D. 
radiodurans, very few references in the scientific literature mention 
the presence of shorter protein fragments after radiation exposure 
in any organism. 

Interestingly, Parodi-Talice and collaborators [34] observed a 
similar pattern in T. cruzi for the proteins glutamate dehydroge- 
nase (GluDH), HSP70, and alpha and beta-tubulins, where lower 
molecular weight isoforms were differentially expressed during 
metacyclogenesis when compared with isoforms with the predicted 



molecular weight. The transformation of epimastigotes into 
metacyclic trypomastigotes is a complex process of diflFerentiation, 
requiring a controlled production of various proteins [34]. 
Similarly, a quantitative time-course proteome analysis for the 
schizont-stage of Plasmodium falciparum (34 to 46 hours after 
invasion) demonstrated that actin-I, enolase, HSPs, and eukaryotic 
initiation factor 4A and 5A presented more than one isoform. The 
isoforms also showed different expression patterns at the different 
time points analyzed. P. falciparum is characterized by a complex 
life cycle, undergoing extensive morphological and metabolic 
changes, which reflects its capacity to survive in different host 
environments [35]. According to the authors, post-translational 
modifications may be a very important strategy for the parasites to 
control gene expression during differentiation [34,35]. 

Differentially Expressed Proteins after Gamma Radiation 
Exposure 

Regarding the differentially expressed proteins, many of the 
listed proteins in Table 2 and Figure S3 are related to the protein 
synthesis process that seems to be upregulated, except for some 
protein spots of the elongation factor 2 that show a reduction in 
their levels. This may be a response to compensate for the 
processing of proteins that occurs after irradiation. This response 
may also enhance the synthesis of specific proteins that wiU 
possibly play a role in the stress response. The results obtained 
from the analyses of translation inhibition and proteomic profile 
after irradiation place de novo protein synthesis as an important 
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Figure 5. Boxplots of peptide molecular weights. A) Distribution of the observed molecular weight in downregulated (green) or upregulated 
(red) protein spots at each time point analyzed. B) comparison between the distribution of the expected (E) and observed (O) molecular weights 
among downregulated or upregulated protein spots 24 and 96 hours after irradiation. A single asterisk corresponds to p<0.05 and a double asterisk 
corresponds to p<0.001. 
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cellular response to gamma radiation. The same pattern is 
observed in D. radiodumns, where proteins related to translation/ 
folding displayed either enhanced or de novo expression in the first 
hour of post-irradiation recovery. Proteins involved with DNA 
repair and oxidative stress alleviation were also induced in D. 
radiodumns under ionizing radiation stress [16]. 

Proteins involved in protein folding processes, such as chaper- 
ones, are mosdy downregulated post-irradiation (Figure S3). This 
represents an unexpected result, since these proteins are classically 
involved with stress response by stabilizing newly synthesized 
protein molecules. Nevertheless, this result is in agreement with 
transcriptomic data observed in microarray experiments [7]. It is 
worth noting that, although HSPs are mosdy downregulated, 
processed forms of these molecules are upregulated and may even 
be functional. Interestingly, the two chaperones localized in the 
endoplasmic reticulum (calreticulin and protein disulfide isomer- 
ase) are upregulated after gamma radiation exposure, which may 
indicate an important role of this compartment in the ionizing 
radiation stress response, suggesting the existence of an unfolded 
protein response-like in this condition [36]. 

Another unexpected result is the downregulation of proteins 
involved in the ATP metabolism (namely the beta subunit of ATP 
synthase and the subunit IV of cytochrome c oxidase), although 



another member of this class is upregulated (cytochrome c oxidase 
subunit V). The outcome of this result is not clear and a more in- 
depth study of the cell energy metabolism would be important. 

Perhaps the most remarkable observation in the post-irradiation 
proteome investigated here is the putative decline in the activity of 
the glycolytic and amino acid metabolism pathways. Several 
important enzymes of glycolysis were downregulated after gamma 
radiation exposure. Accordingly, the only enzyme (pyruvate 
phosphate dikinase) from gluconeogenesis hsted here was upregu- 
lated. Most enzymes involved in the amino acid metabolism were 
also downregulated, but shorter isoforms of the GluDH were 
upregulated after irradiation. They consist of three isoforms with 
experimental molecular weights (15 kDa) lower than the predicted 
values (45 kDa), suggesting once again the occurrence of post- 
transcriptional modifications/processing of important metabolic 
enzymes during the stress response. GluDH catalyzes the NAD- 
and/ or NADP-dependent reversible deamination of L-glutamate 
to form alpha-ketoglutarate and is essential for the metabolism of 
amino nitrogen in organisms ranging from bacteria to mammals 
[37]. T. cruzi has a metabolism that is largely based on die 
consumption of amino acids, mainly, proline, aspartate, and 
glutamate, which constitute the main carbon and energy sources of 
the epimastigote forms. In T. cruzi, GluDH has NADP-specilic 
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activity [38], indicating that it may serve as a pentose-phosphate 
shunt-independent source of NADPH in these parasites. Taken 
together, these results suggest that the parasite experiences an 
overall reduction on its energy metabolism as a consequence of its 
growth arrest after irradiation. 

We have identified four proteins classified as redox sensors in 
this study. While two of these are downregulated (both oxidore- 
ductases), the other two are upregulated and these are both 
tryparedoxins, which efficiendy reduce hydrogen peroxide [39]. 
Throughout its life cycle, T. cruzi is exposed to various stresses in 
different environments: the invertebrate (triatomine bugs) and the 
vertebrate hosts. One of the most deleterious consequences of 
oxidative stress may be the formation of DNA lesions. Guanine is 
the most susceptible base to oxidation, due to its low redox 
potential, and the 7,8-dihydro-8-oxoguanine (8-oxoG) is the most 
common lesion. When 8-oxoG is inserted during DNA replication, 
it can generate double-strand breaks, which makes this lesion 
severely deleterious. Recently Aguiar et ai, 2013, demonstrated 
that parasites overexpressin MutT are more resistant to the 
oxidative stress caused by hydrogen peroxide (H2O2) treatment. 
The MutT enzyme product, 8-oxod-GMP, can generate an 
oxidative stress signal, enabling the cells to overcome this stress. 
MutT hydrolyses 8-oxo-dGTP in the nucleotide pool, returning it 
to the monophosphate form so that it cannot be incorporated into 
DNA by polymerases. Parasites overexpressing heterologous 
MutT also increase the levels of cytosolic and mitochondrial 
peroxidases (TcCpx and TcMPx) after H2O2 treatment. Taking 
this into account and also that parasites subject to gamma 
radiation experience oxidative stress and increase the levels of 
some antioxidant enzymes not immediately after irradiation, but 
later after irradiation, we could suggest that T. cruzi does not 
respond directly to ROS production as a consequence of 
irradiation, but to 8-oxo-dGMP that is generated subsequently. 
The nucleotide 8-oxo-dGMP, or another secondary metabolite 
generated from this process, could be acting as a second messenger 
to the cell and indicating the presence of oxidative stress. 

Recendy, Krisko & Radman proposed a new paradigm when a 
cell is subject to ionizing radiation: the proteome rather than the 
genome is the primary target in radiation-induced cell death. 
This paradigm has been supported by several experimental 
evaluations showing that D. radiodurans has a way of protecting its 
proteins from oxidative damage [40] . Indeed, a strong correlation 
between intracellular Mn/Fe concentration ratios and bacterial 
resistance to radiation has been shown, in which the most 
resistant bacteria tolerates 300 times more Mn^""" and three times 
less Fe^^ than the most radiation-sensitive bacteria [1.^]. 
Manganese ions prevent the formation of iron-depend(uit ROS 
through the Fenton reaction, acting as chemical antioxidant 
protectors. Furthermore, measurements of protein carbonyl 
groups in D. radiodurans revealed that Mn^"*" accumulation 
prevented protein oxidation; these results were also observed in 
other radioresistant bacteria [41-44]. 

Furthermore, the level of oxidative protein damage caused 
during irradiation controls the survival of many organisms [Bdelhid 
rotifers, a class of freshwater invertebrates, Caenorhabditis elegans, the 
bacteria D. radiodurans, and Halobacterium salinarum), which are 
extremely resistant to ionizing radiation [44-47]. The principal 
factor responsible for this extraordinary radioresistance is their 
great antioxidant protection of their cellular constituents, includ- 
ing those required for DSB repair, allowing them to recover from 
stress and continue reproduction [46]. 

An important finding of this study is the significant upregulation 
of three hypothetical proteins after gamma radiation. This may 
indicate a role for species-specific proteins in the response to stress 



after ionizing radiation exposure, since these most likely represent 
proteins with no homologues in other species. Similarly, in an 
initial D. radiodurans proteome study, hypothetical proteins were 
identified and further proved to be crucial for the response to 
radiation in this bacterial species [47]. 

Sghaier and collaborators have recently published a study on 
the amino acid composition of proteins from radiation-resistant 
bacteria [48]. The authors report that such proteins bear more 
small amino acids and fewer aromatic rings. We have also assessed 
the amino acid composition of T. cruzi proteins in a slightiy 
dilferent perspective. Amino acid counts were performed for 
upregulated T. cruzi proteins after gamma radiation exposure (we 
have considered as upregulated the proteins that were more 
abundant than in NI cells at least in one time point) and for the 
orthologues in T. brucei of T. cruzi upregulated proteins. In both 
cases, amino acid counts were normalized by the count performed 
in the set of all annotated proteins of the respective Trypanosoma sp. 
The hypothesis was that proteins with important roles after irradia- 
tion in T. cruzi would have an amino acid composition different than 
what is observed in the set of all T. cruzi proteins and in the 
respective orthologues in T. brucei (which is not radio-resistant). 

When we compared T. cruzi proteins that were upregulated 
after radiation exposure with the entire set of annotated proteins 
from this parasite, we observed that the former have in general 
fewer polar, hydrophobic, and small amino acids (although some 
amino acids in these classes are more frerjuent). In addition, 
upregulated proteins have fewer aromatic amino acids (except for 
tyrosine, which is more frequent) and less sulfur-containing cysteine 
residues. 

Conclusions 

Using 2D-DIGE and MS, we have identified 543 protein spots 
differentially express(-d after gamma radiation exposure. The 
presence of multiple isoforms was observed for more than half of 
the identified proteins, most of which are shorter than the 
annotated protein size in the T. cruzi genome. Additionally, there 
was a strong correlation for lower molecular weight peptide spots 
to be overexpressed. This result could be explained by de novo 
protein synthesis of different isoforms, protein processing, and/ or 
modification events subsequent to radiation exposure. This 
observation indicates that post-translational control of gene 
expression have an important role in the parasite response to 
gamma radiation stress. The inhibition of protein synthesis in face 
of gamma radiation was shown to have a significant effect 
decreasing parasite growth and survival rates, highlighting the 
importance of active translation for parasite recovery after 
exposure to ionizing radiation. 

We have annotated all 53 proteins identified by MS according 
to their biological roles. Several proteins were represented by 
multiple spots, and most of them had molecular weights lower 
than predicted. As a consequence of this observation, we cannot 
precisely state which biological processes are upregulated versus 
downregulated, since the different protein isoforms may not 
function in the same way as the full-length protein. Nevertheless, 
some tendencies could be observed in this study, including 
changes in the following biological processes: upregulation of the 
protein synthesis process, downregulation of protein folding 
(except for the upregulation of two endoplasmic reticulum 
chaperones), downregulation of the ATP generation j)athway, 
glycolysis, and amino acid metabolism, and the upregulation of 
two tryparedoxins (which reduce hydrogen peroxide in response 
to oxidative stress). 
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Finally, taking into account the translation inhibition results 
obtained herein, together with the observed proteomic profile after 
irradiation, we can conclude that de novo protein synthesis is an 
essential cellular response to gamma radiation. 

Supporting Information 

Figure SI Electrophoretic analysis of total protein extracts of 
irradiated and NI epimastigote cells. Total protein extracts were 
obtained for each time point NI, 4, 24, and 96 hours after 
irradiation. Samples were subjected to 12% SDS-PAGE and 
stained with coomassie blue. 
(TIF) 

Figure S2 Differentially expressed isoforms of prostaglandin F2 
alpha synthase. The upregulated protein spot (161) shows a lower 
molecular weight when compared with the downregulated 
proteins spots (14, 111, 113, 114, and 144). 
(TIF) 
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